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HE 要 : 为 定量 区 分 退耕 还 林 还 草 背 景 下 北方 典型 农 牧 交错 带 植 被 变化 过 程 中 气候 变化 和 人 类 活 


动 的 相对 贡献 率 , 以 宁夏 盐池 上 


2000—2020 年 植被 变化 为 研究 对 象 ,基于 MODIS13Q1-NDVI 数 


据 、 地 表 和 覆盖 数据 和 气象 数据 ,采用 Thornthwaite 纪念 模型 和 CASA (Carnegie-ames-stanford ap- 
proach ) 模 型 分 别 估算 了 逐年 的 潜在 净 初 级 生产 力 (Potential net primary productivity, PNPP) 和 实际 
净 初 级 生产 力 (Actual net primary productivity, ANPP ) 。 综 合 运用 趋势 分 析 、 相 关 分 析 和 差 值 比较 法 
分 析 了 2000 一 2020 年 盐池 县 净 初 级 生产 力 (Net primary productivity, NPP) 时 空 变化 特征 及 其 驱动 
A ,并 定量 确定 了 气候 因子 和 人 类 活动 对 盐池 县 植被 变化 的 相对 贡献 率 。 结 果 表明 :(1) 在 2000 一 


2020 年 盐池 县 NPP 总 体 旦 上 升 趋势 ,但 存在 着 显著 的 空间 异 质 性 ,主要 表现 为 植被 NPP dr dE IX 
R, 上 且 改 善 或 退化 程度 也 存在 显著 的 空间 分 异 。 植 被 改善 区 主要 分 布 于 荒 
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地 带 ,(2) 植被 变化 归 因 分 析 表 明 , 在 植被 改善 


区 ,气候 变化 和 人 类 活动 


# 同 主导 驱动 了 植被 的 改 


善 , 但 气候 变化 的 相对 贡献 率 ($9.77% ) 大 于 人 类 活动 的 相对 贡献 率 (40.23% ) , 而 在 植被 退化 区 ,人 
类 活动 的 相对 贡献 率 (91.77%) 则 显著 高 于 气候 变化 的 相对 贡献 率 (8.23%)。(3) 驱动 力 分 析 表 明 ， 
研究 区 植被 NPP 变 化 与 降水 量 呈 显著 正 相 关 ,而 与 气温 的 相关 性 较 弱 ;而 人 类 活动 是 驱动 植被 退 
化 区 NPP 下 降 的 主要 原因 ,但 负 向 影响 力 有 所 减弱 。 总 体 而 言 ,气候 变化 是 植被 改善 区 的 主要 驱动 


力 , 而 人 类 活动 是 植被 退化 区 的 主要 驱动 因素 , 


表 者 共同 作用 则 使 盐池 县 整体 生态 环境 得 以 改善 。 


X 键 词 : 净 初 级 生产 力 ; 荒漠 草原 ; 植被 变化 ; 驱动 力 ; 盐池 县 
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盐池 县 地 处 我 国 北方 农 牧 交错 带 ,在 近 几 十 年 
内 ,在 全 球 气候 变化 背景 下 ,由 于 过 度 放 牧 、 倒 山 种 
地 矿产 开采 药材 (甘草 ) 采 集 、 编 制 原料 的 樵 采 等 人 
类 活动 导致 植被 退化 .土地 沙化 和 水 土 流失 加 重 "。 
为 逆转 沙漠 化 进程 .修复 退化 草原 ,盐池 县 于 2002 
年 启动 退耕 还 林 项 目 , 于 2003 年 开始 推行 草原 禁 牧 
封 育 政策 和 退 牧 还 草 工 程 ,于 在 控制 土壤 侵蚀 Be 
高 生态 环境 质量 恢复 草原 植被 .提高 植被 覆盖 度 
和 遏制 土地 沙漠化 进程 。 同 时 ,气候 格局 的 演变 促 


的 相对 贡献 ,对 深入 理解 植被 变化 驱动 机 制 、 保 障 
草地 生态 安全 和 控制 荒漠 化 发 展 具有 重要 意义 。 
因此 , 亚 需 构建 植被 变化 动态 监测 方法 并 辨析 各 影 
响 因 子 的 相对 贡献 率 。 

植被 净 初 级 生产 力 (Net primary productivity, 
NPP) 是 植物 通过 光合 作用 转换 累积 的 有 机 物质 总 
量 , 是 表征 植物 活动 的 重要 变量 ,可 以 准确 反映 植 
被 生长 状况 ”。 由 于 NPP 对 气候 变化 和 人 类 活动 敏 
感 ,是 监测 草地 变化 、 区 分 气候 变化 和 人 类 活动 影 


进 了 生态 工程 成 效 的 提升 和 植被 的 恢复 ”。 气 候 变 
化 和 人 类 活动 是 植被 变化 的 两 大 主要 驱动 因子 ” ， 
如 何 定量 区 分 气候 变化 与 人 类 活动 在 植被 变化 中 
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响 的 可 靠 指标 “”"。 目 前 定量 区 分 气候 变化 与 人 类 
活动 在 植被 变化 中 相对 贡献 率 的 主要 方法 有 回归 
模型 法 、 差 值 比较 法 、 基 于 生物 物理 过 程 的 模型 方 
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法 ”。 其 中 ,将 差 值 比较 法 与 基于 NPP 的 评价 方法 
相 结 合 已 成 为 植被 变化 中 将 气候 变化 因素 和 人 类 
活动 因素 区 分 研究 中 应 用 最 多 的 方法 ”“。 其 原理 如 
下 :假设 未 受 人 类 活动 干扰 的 、 仅 以 气候 为 主要 驱动 
力 的 植被 NPP 为 潜在 净 初 级 生产 力 (Potential net 
primary productivity, PNPP)” ,主要 借助 气候 模型 估 
算 丰 ;实际 净 初 级 生产 力 (Actual net primary produc- 
tivity, ANPP) 则 为 气候 变化 和 人 类 活动 共同 影响 下 
的 植被 NPP, 可 借助 CASA ( Carnegie- ames- stanford 
approach ) 模 型 进行 估算 , 因 其 所 需 参 数 少 旦 易于 获 
取 ” ,是 目前 应 用 最 广泛 的 估算 模型 ;而 PNPP 与 
ANPP 的 差 值 被 认为 是 人 类 活动 作用 下 的 NPP 损 失 
或 增加 量 , 即 人 类 活动 影响 下 NPP 的 变化 值 (Hu- 
man net primary productivity, HNPP) John 等 成功 
应 用 差 值 比较 法 区 分 了 人 类 活动 和 降水 对 蒙古 高 
原 植 被 覆盖 变化 的 影响 ,结果 表明 植被 NPP 下 降 的 
原因 主要 是 源 自 放牧 强度 和 人 口 密度 的 增加 ,而 植 
被 的 上 升 则 归 因 为 生态 恢复 工程 及 围 封禁 牧 政 策 
的 实施 。 田 海 静 "和 Piao 等 采用 差 值 比较 法 定 
量 评估 了 我 国人 为 生态 恢复 工程 的 成 效 ,认为 除 气 
候 变 化 外 ,林业 建设 工程 对 植被 变化 有 显著 贡献 。 
说 明 差 值 比 较 法 不 仅 可 以 厘定 各 类 影响 因子 贡献 
率 的 区 域 差异 ,也 可 以 量化 评价 人 类 活动 在 NPP 变 
化 过 程 中 的 相对 贡献 率 。 这 相 较 于 以 往 气候 变化 
和 人 类 活动 整体 性 分 析 “” ,能 更 好 地 解释 区 域 植 
被 变化 的 空间 分 异 及 其 驱动 机 制 。 因 此 ,本 研究 充 
分 考虑 了 气候 变化 和 人 类 活动 因素 的 空间 蜡 质 性 
以 及 气象 数据 的 可 获取 性 ,在 县 域 尺 度 上 探讨 宁夏 
典型 荒漠 草原 分 布 区 在 大 型 生态 工程 (退耕 还 林 还 
草 工程 ) 和 围 封 禁 牧 政策 实施 前 后 的 植被 变化 及 其 
驱动 机 制 ,以 植被 NPP 为 量化 指标 ,利用 CASA 模 型 
估算 2000 一 2020 年 盐池 县 ANPP 的 动态 变化 ,并 结 
合 基于 气候 模型 估算 的 PNPP 及 二 者 差 值 ,定量 评 
佑 气候 变化 与 人 类 活动 在 植被 动态 变化 过 程 中 的 
相对 贡献 率 及 其 时 空 变化 。 


1 研究 区 概况 


盐池 县 位 于 宁夏 东部 (37*04'~38?10'N ,106*30'~ 
107°47'E) ,土地 面积 6749 km’, 占 宁夏 总 面积 的 
10.16%。 地 形 上 ,属于 鄂尔多斯 台地 向 黄土 高 原 的 
过 渡 带 ;气候 上 ,属于 干旱 向 半 干 旱 区 及 草原 向 藻 
漠 草 原 区 的 过 渡 带 ,具有 典型 的 中 温带 大 陆 性 气候 


特征 ,常年 干旱 少雨 , 风 大 沙 多 ;资源 分 布 上 , 则 是 
宁夏 矿产 资源 的 集中 分 布 区 之 一 。 研 究 区 年 降水 
量 250-350 mm, 年 均 气 温 8.1 7C ,年 蒸发 量 高 达 
2200 mm ,是 年 降水 量 的 8.8 倍 左右 , 且 降 水 量 多 集 
中 在 7 一 9 月 ,盛行 风 为 西风 和 西北 风 。 地 带 性 土壤 
以 灰 钙 土 ` 风沙 土 为 主 ,结构 松散 .肥力 较 低 , 易 受 风 
蚀 和 水 蚀 的 影响 。 参 照 1:1000000 中 国 植被 图 集 
可 将 当地 主要 植被 类 型 大 致 分 为 6 类 ,分 别 为 典型 
草原 .荒漠 草原 .荒漠 、 粮 食 及 经 济 作物 落叶 灌 从 
和 草 多 (网 1)。 


mu 粮食 及 经 济 作物 
mm 落叶 灌 丛 
mu Ei 


图 1 研究 区 植被 类 型 
Fig. 1 Vegetation types in the study area 


2 数据 与 方法 


2.1 数据 来 源 

2.1.1 i£ E d& d& — 本 研究 采用 归 一 化 植被 指数 
(Normalized difference vegetation index, NDVI) 和 
NPP 2 个 反映 植被 生长 状态 的 指标 数据 来 分 析 人 研究 
区 的 草地 植被 时 空 动态 , 均 来 自 于 Level-1 and At- 
mosphere Archive & Distribution System Distributed 
Active Archive Center(LAADS DAAC) jy, rp, 
NDVI 数 据 集 选用 时 间 序 列 为 2000 一 2020 年 的 MO- 
DIS 陆地 标准 产品 MOD13Q1 数 据 集 (空间 分 辨 率 为 
250 m; 时 间 分 辩 率 为 16d), 使 用 最 大 值 合成 法 (Max- 
imum value composition , MVC) 生 成 逐 月 最 大 NDVI 数 
Ji" , Æ CASA 模型 估算 ANPP 的 重要 参数 ”。 同 
时 ,选用 同时 间 序 列 的 MODIS 陆地 标准 产品 
MOD17A3HGF 数据 集 (空间 分 状 率 为 500 m; 时 间 分 
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辩 率 为 1a) 用 于 验证 CASA 模 型 NPP 模 拟 精 度 。 
2.1.2 气象 数据 气象 数据 来 自 于 国家 气象 信息 中 
心 中 国 气 象 数 据 网 (http://data.cma.cn/) ,本 研究 
选取 盐池 及 其 附近 14 个 气象 站 点 数据 ,包括 :月 降 
水 量 月 平均 气温 日照 百分率 及 经 纬度 信息 ,时间 
跨度 为 2000 一 2020 年 。 根 据 气象 站 点 经 纬度 信息 ， 
基于 ArcGIS 10.2 的 Geostatistical Analyst 模块 中 的 
克 里 金 插值 法 对 气象 数据 进行 插值 运算 。 
2.1.3 4, 7E E dk de RAMA APRA 
高 技术 研究 发 展 计划 (863 计划 ) 全 球 地 表 履 盖 遥 感 
制图 与 关键 技术 研究 项 目的 重要 成 果 GlobeLand30 
数据 集 ,空间 分 辨 率 为 30 m, 时 间 序 列 为 2000、 
2010.2020 4. 
模型 中 所 有 输入 参数 , 均 转 换 为 Albers equal 
area conica- WGS 84 投 影 坐 标 系统 , 重 采 样 为 与 ND- 
VI 数 据 相 同 的 空间 分 辩 率 , 且 裁 剪 至 像 元 行列 数 
相同 。 
2.2 指标 计算 方法 
2.2.1 植被 实际 净 初 级 生产 力 (ANPP) 本 人 研究 采 
用 朱文 泉 等 2 的 CASA 改进 模型 对 盐池 县 2000 一 
2020 年 的 植被 ANPP 进行 估算 ,该 模型 可 对 生态 系 
统 的 生产 力 .土壤 碳 进行 模拟 ,并 已 广泛 应 用 于 各 
RET NPP fi rp, CASA 模型 为 NDVI 地 表 
覆盖 类 型 .气象 数据 共同 驱动 的 光 能 利用 率 模 型 ， 
可 通过 植被 吸收 的 光合 有 效 辐射 (APAR) 和光 能 利 
用 率 (g) 来 确定 , 即 : 
ANPP(x, 2) = APAR(x,t) X &(x,t) (1) 
IEP ANPP(x,t), APAR(x,) .e(x,t)2 9] AJ 2C x TE t A 
的 植被 初级 生产 力 (g Cm?) 吸收 的 光合 有 效 辐射 
(MJ* m^) RISE BR CBE HAS (g C MI); 
APAR(x,t)= SOL(,) xX FPAR(x,1)X0.5 — (2) 
IP : FPAR(x,t) 4 FE BERT 2G 6r A 23058 9] 0] HR AC EG 
例 ;常数 0.5 为 植被 利用 0.38~0.71 nm 波长 的 太阳 有 
效 辐射 占 太 阳 总 辐射 比例 ;SOLC) 为 像 元 * 在 :月 
的 太阳 总 辐射 量 (MJ.m2?) ,参照 和 清华 等 ”提出 的 
中 国 西部 地 区 太阳 总 辐射 计算 公式 计算 : 
SOL(x,)- SOL,(,) x(0.185 +0.595S) (3) 
式 中 :SOLo(w,) 为 像 元 x 在 1 月 的 天 文 辐射 (MJ.m”);S 
为 日 照 百分率 。 
FPAR 人 取决 于 植被 类 型 和 植被 的 覆盖 状况 ， 
与 NDVI 具 有 良好 的 线性 关系 , 即 : 


FPAR(x,t)= 
[NDVI(x,) - NDVI. X(FPAR,,, — FPAR,,,) 
NDVI, - NDVI, 


*FPAR,.. 


(4) 
IEP NDVI (x,t) NAIC x YE t H B9 NDVI ME ; NDV Linas 
NDVL,。 分 别 为 植被 NDVI 最 大 值 和 最 小 值 ;FPAR,、 
FPAR, 的 取 值 与 植被 类 型 无 关 , 分 别 为 0.95 .0.001。 
elx, t) =T (xt X T o(x.t) Xx W, (x,t) X E nax (5) 
WP ex.) NR x E LH SEE RE Al HK Cg C- 
MJ) ; T. (x.t) I T. xut) Jg COSE DEBERI AS H ; 
W.(x,t) 97K 4 GE S Us] ZR BL 5 Ena AT PRAEDAS. BY Se 
大 光 能 利用 率 ,sw 的 取 值 参 考 Zhu 等 2 提出 的 各 植 
被 类 型 区 的 最 大 光 能 利用 率 模拟 值 。 
2.2.2 植被 潜在 净 初 级 生产 力 (PNPP) 本 研究 应 
用 Thornthwaite 纪念 模型 中 进行 植被 PNPP 的 估 
算 。 该 模型 是 被 广泛 接受 的 PNPP 估算 模型 ,是 基 
于 实地 测量 的 NPP 数 据 与 研究 区 气温 和 降水 量 数 
据 之 间 的 最 小 二 乘 回归 得 出 的 。 计 算 公式 如 下 


PNPP = 3000 x [l 2 p M (6) 

J= 1.05r - (7) 
1+(1+1.05r/L) 

L= 3000 + 25t + 0.05¢° (8) 


式 中 :PNPP 为 年 总 潜在 NPP(g Cm?) o WEF 
际 蒸 散发 量 (mmy) ;7 为 年 总 降水 量 (mm) ;为 年 平 
JR BU E (mm) ;为 年 平均 气温 (% )。 
2.2.3 差 值 比较 法 ”基于 NPP 的 差 值 比 较 法 能 直 
观 、 真 实地 反映 生态 环境 中 植被 的 变化 情况 ,从 而 
确定 各 自 的 相对 作用 ,目前 这 一 方法 是 研究 植被 变 
化 驱动 机 制 的 新 途径 ”% , HNPP Jy PNPP 5j ANPP 的 
差 值 , 即 : 

HNPP = PNPP - ANPP (9) 
式 中 :PNPP 为 潜在 净 初 级 生产 力 ;ANPP 为 实际 净 
初级 生产 力 ;HNPP 为 人 类 活动 影响 下 NPP 的 变化 
值 ,HNPP 值 的 正 负 为 人 类 活动 对 NPP 影 响 的 方向 ， 
正 值 表示 人 类 活动 影响 下 NPP 减少 , 负 值 表示 人 为 
影响 下 NPP 增 加。 为 辩 析 人 类 活动 对 植被 NPP 的 
影响 趋势 和 方向 ,本 文 参照 周 妍 妍 等 的 研究 方 
法 ,将 HNPP 的 正 负 与 变化 趋势 进行 全 加 分 析 , 从 而 
获取 人 类 活动 对 植被 NPP 的 影响 趋势 及 其 空间 分 布 。 
2.2.4 NPP 年 际 变化 率 ”NPP 动态 可 直接 反映 植被 
的 改善 或 退化 状况 ,ANPP 的 变化 趋势 可 以 监测 植 
被 变化 ,PNPP 和 HNPP 的 变化 趋势 可 以 评估 气候 变 
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化 和 人 类 活动 对 植被 变化 的 相对 影响 。 本 研究 采 
用 基于 最 小 二 乘法 的 一 元 线性 回归 分 析 方 法 计算 
2000 一 2020 年 3 种 定义 下 的 NPP 年 际 变化 趋势 。 公 
式 如 下 : 


nx YExNPP)- yiyNPP 
i= i=l d 


nx » 2:2 
式 中 :n 为 研究 时 间 段 内 的 总 年 份 数 ; NPP 为 第 i 年 
的 NPP 值 ;Slope 为 NPP 的 变化 趋势 ,计算 结果 为 正 ， 
NPP 有 增加 趋势 。 采 用 下 检验 验证 变化 趋势 的 显著 
性 ,检验 结果 仅 表示 变化 趋势 的 可 置信 程度 ,与 变 
化 速率 无 关 , 公 式 如 下 : 


Slope = (10) 


(£, - à 
F=(n-2) — (11) 


2-5 =a) 


式 中 :为 第 i 年 NPP 的 回归 值 ; x 24 21 ala] NPP Æ 
均值 ; x, 为 第 ;年 的 NPP 值 ,z=21。 基 于 下 检验 可 以 
将 结果 分 为 以 下 6 种 等 级 : 极 显 著 改 善 (Slope>0,P< 
0.01) .显著 改善 (Slope>0,0.01< Ps 0.05) .无 显著 改 
善 (Slope>0, P>0.05) .无 显著 退化 (Slope<0, P» 
0.05) .显著 退化 (Slope<0,0.01< P< 0.05) . 极 显 著 退 
化 (Slope<0,P<0.01)。 

人 研究 时 间 段 内 NPP 的 变化 量 ( ANPP ) ,用 以 下 
AG : 


ANPP =(n - I) x Slope (12) 
2.2.5 NPP 变 化 娄 因 定量 评估 方法 ”为 了 对 气候 变 
化 和 人 类 活动 在 植被 的 改善 或 退化 中 相对 贡献 率 
进行 区 分 ,可 以 通过 计算 3 种 定义 下 的 NPP 年 际 变 
化 趋势 Saner, Sonera Suv» LA Be NPP 变化 量 APNPP、 
AHNPP 来 实现 ,根据 刘斌 等 相关 NPP 变 化 相对 贡 
献 的 情景 设 定 方案 ,定义 以 下 6 种 情景 ( 表 1)。 参 照 


表 1 不 同情 景 下 气候 变化 和 人 类 活动 在 植被 变化 中 相对 贡献 的 评估 方法 


Tab. 1 Methods for assessing the relative contributions of climate change and human activities 


to vegetation change under different scenarios 


植被 变化 状况 情景。 Snw Same ”气候 变化 的 相对 贡献 /96 人 类 活动 的 相对 贡献 /%6 说 明 
植被 改善 。 ”情景 ! >0 >0 100 0 气候 变化 主导 植被 改善 
(Sww>0) 情景 <0 <0 0 100 人 类 活动 主导 植被 改善 
情景 3 >0 «0 IAPNPP| 100g, . . lAHNPP!  zji10ge, 共同 主导 植被 改善 ,以 其 各 自 变化 
IAPNPPI+IAHNPPI IAPNPPI + IAHNPPI 量 所 占 的 比例 为 各 自 的 相对 贡献 
植被 退化 ”情景 4 <0 — «0 100 0 气候 变化 主导 植被 退化 
(Swm<0) 情景 5 >0 >0 0 100 人 类 活动 主导 植被 退化 
情景 6 <0 »0 IAPNPP| —— x100g IAHNPP| ggg, 共同 主导 植被 退化 ,以 其 各 自 变化 


IAPNPPI + IAHNPPI 


表 1 中 各 类 情景 ,可 分 别 计 算出 人 类 活动 .气候 变化 
对 植被 变化 的 相对 贡献 率 。 


3 结果 与 分 析 


3.1 ANPP 的 模拟 及 其 验证 

本 研究 从 像 元 尺度 对 CASA 模型 估算 2000 一 
2020 年 ANPP 模 拟 值 和 同时 段 的 MODIS NPP 产 品 
进行 一 一 对 应 精度 验证 ”与 相关 性 分 析 ( 图 2)。 结 
果 表 明 ,ANPP 估 算 结果 和 MODIS NPP 值 基本 吻合 ， 
模拟 精度 较 高 ,决定 系数 达 0.82(P<0.01) ,上 且 估算 结 
果 比 MODIS NPP 值 略 低 。 
3.2 ANPP 的 时 空 变化 趋势 

由 图 3 可 以 看 出 ,2000 一 2020 年 盐池 县 年 均 


IAPNPPI* [AHNPPI 量 所 占 的 比例 为 各 自 的 相对 贡献 


UE Sue Ser Su J ANPP PNPP , HNPP 的 年 际 变 化 趋势 ;APNPP .AHNPP 为 PNPP .HNPP 的 变化 量 。 


600 r 
y=1.49x-19.43 
p R?=0.82 
ai P«0.01 
g 
o 400r 
= 
x 
r3 200 上 
2 
0 L 1 L 


f 
100 200 300 400 
MODIS NPP 值 /g C-m?-a* 


注 :ANPP 为 实际 净 初 级 生产 力 ;NPP 为 净 初 级 生产 力 。 下 同 。 
图 2 MODIS NPP 和 ANPP 模 拟 值 对 比 
Fig. 2 Comparison of MODIS NPP and 


ANPP simulation values 
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tidwi 

. EM ANPP 密度 整体 上 呈现 出 上 升 趋势 ,增长 速率 为 
ixl . 5.59 g C+ m?- a`! (770.54, P«0.01) ,表明 盐池 县 
E 200 | mo 2 ANPP 总 体 上 呈 增 加 趋势 。 但 从 不 同时 间 段 来 看 ， 
È iL. [A y-5.59:140,53 整体 上 呈现 波动 中 上 升 趋势 。 按 照 植被 波动 方向 
< 100. ror 来 看 ,其 中 在 2004—2005 ,2012—2013 ,2018—2020 

50 — ——— — — — — 年 出 现 小 幅 的 负 向 波动 。 
FS S ESSE MSS 2000 一 2020 年 盐池 县 年 均 ANPP 空间 分 布 呈 从 
Eh 东南 回 西 北 逐 渐 减 少 的 趋势 ,空间 上 存在 几 个 明显 
图 3 2000_2020 年 盐池 县 植被 ANPP 的 年 际 变化 no a dn 
Fig.3 Interannual ANPP change in Yanchi 要 分 布 在 耕地 集聚 区 和 块 状 扬 黄 灌 区 , ANPP 4} i 
County during 2000—2020 250-500 g Cm *…a 之 间 。 随 着 退耕 还 林 等 生态 工 

(a) 2000 一 2020 年 多 年 平均 ANPP (b) 2000 年 ANPP 


NORA 
EDS 
SUE ERS 
Ie? wipe 
D sig MNS 
(c) 201042 ANPP (d) 20204F ANPP 


ibang 


Garris 
Bi 


SiC Fe S 


fee Weer A 


0 10km 0 10km 
LARA [乡镇 边界 «120120 150170190 210230250 280 310 350 400 >400 
ANPP/g C:m?:2 


图 4 盐池 县 多 年 平均 ANPP 及 2000.2010.2020 年 ANPP 的 空间 分 布 
Fig.4 Spatial distributions ofANPP in 2000, 2010, 2020 and multi-year average in Yanchi County 
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程 实施 范围 的 增 大 ,出现 高 值 区 范围 扩大 的 现象 ; 
低 值 区 主要 分 布 在 高 沙 窝 镇 、 大 水 坑 镇 、 麻 黄山 乡 
的 沙 地 区 域 ,ANPP 介 于 30~120 g Cm*…a! 之 间 。 
2000—2020 年 植被 NPP 变化 趋势 的 显著 性 分 
析 表 明 ,盐池 县 植被 NPP 整体 上 呈 显 著 改 善 趋势 
(图 5a), 且 是 改善 趋势 的 面积 比例 (99.39%) 远 大 于 
呈 退 化 趋势 的 面积 比例 (0.61% )。 尽 管 ANPP 总 体 
呈 改 善 趋势 ,但 其 改善 程度 存在 显著 的 空间 分 异 。 
其 中 ,ANPP 呈 极 显著 改善 的 区 域 占 研究 区 面积 的 
68.91% , 呈 显 著 改 善 的 面积 占 17.23%。 就 空间 分 布 
而 言 , 极 显著 改善 的 区 域 遍 布 于 研究 区 全 境 , 呈 退 
化 趋势 的 区 域 则 零星 分 布 于 荒漠 和 荒漠 草原 及 人 
口 密集 区 (图 5b)。 
3.3 NPP 变化 归 因 分 析 
3.3.1 和 气候 变化 和 人 类 活动 对 NPP 改 善 区 的 相对 贡 
ae 从 图 6a 可 以 看 出 :在 植被 改善 区 ,气候 变化 
人 类 活动 共同 主导 作用 的 相对 贡献 率 为 74.14%， 
而 气候 变化 和 人 类 活动 的 相对 贡献 率 分 别 为 
22.65%、3.21%。 从 空间 分 布 上 来 看 ,气候 变化 主导 
驱动 ( 即 相 对 贡献 率 大 于 50% ) 的 植被 NPP 改 善 区 
主要 分 布 于 盐池 县 西部 的 高 沙 窜 镇 、 王 乐 井 乡 的 人 
口 稀 跑 分 布 区 ;人 类 活动 主导 驱动 ( 即 相 对 贡献 率 
大 于 50%) 的 植被 NPP 改善 区 主要 分 布 于 盐池 县 东 


(a) ANPP 年 际 变化 趋势 


dizit 


TIP 


0 10km 


部 的 花 马 池 镇 、 高 沙 窜 镇 青 
密集 区 (图 7)。 
3.3.2 气候 变化 和 人 类 活动 对 NPP 退 化 区 的 相对 贡 
献 率 如 图 6b 所 示 , 人 类 活动 的 相对 贡献 率 为 
72.48% ,而 气候 变化 的 相对 贡献 率 仅 为 0.47% , 而 二 
者 共同 主导 的 相对 贡献 率 为 27.05%。 说 明 人 类 活 
动 是 局 地 植被 退化 的 主要 原因 。 从 空间 分 布 上 来 
看 ,气候 变化 主导 植被 退化 的 区 域 面积 占 比 仅 为 
8.23% ,而 91.77% 的 植被 退化 区 均 为 人 类 活动 引起 
主要 分 布 于 汉 记 沟 乡 (图 8)。 
3.3.3 和 气候 变化 和 人 类 活动 对 NPP 影响 程度 宝 间 分 
异 及 其 相对 贡献 率 ， 在 植被 改善 /退化 区 2 个 因素 
的 相对 贡献 率 表 现 出 明显 的 空间 异 质 性 。 在 植被 
极 显 著 改 善 区 ,人 类 活动 对 植被 改善 的 相对 贡献 率 
H 76.0196 ,气候 变 化 的 相对 贡献 率 为 23.99% ;在 植 
被 显著 改善 区 ,气候 变化 对 植被 改善 的 相对 贡献 率 
dies 而 人 类 活动 的 相对 贡献 率 为 17.87%( 表 
)。 从 不 同 植被 类 型 上 看 ,人 类 活动 为 荒 江 和 匾 漠 
,气候 变化 为 典型 
草原 .粮食 及 经 济 作物 . 草 多 .落叶 灌 丛 等 区 域 植被 
改善 的 主要 因素 ( 表 3)。 在 各 等 级 植被 退化 区 ,人 
类 活动 对 植被 退化 的 相对 贡献 率 均 达 到 70% 以 上 
( 表 2)。 从 不 同 植被 类 型 上 看 ,植被 退化 区 仅 有 荒 


乡 . 王 乐 井 乡 的 人 口 


(b) 变化 趋势 显著 性 分 析 


行政 边界 
O 县 界 
O 乡镇 边界 

显著 性 检验 结果 


图 5 2000 一 2020 年 盐池 县 ANPP 年 际 变化 趋势 及 显著 性 检验 
Fig. 3 Interannual variation trend and significance test of ANPP in Yanchi County during 2000—2020 
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(a) 植被 改善 区 


"d 


行政 边界 
co BF 

C3 乡镇 边界 
情景 空间 分 布 

mu 人 类 活动 主导 
mm 气候 变化 主导 
mm 共同 主导 


(b) 植被 退化 区 


SZ 


va 高 沙 窜 镇 


图 例 
行政 边界 
mE» 

L3 乡镇 边界 
情景 空间 分 布 

m 人 类 活动 主导 
Em 气候 变化 主导 
mm 共同 主导 


图 6 气候 变化 .人 类 活动 及 二 者 共同 作用 在 植被 改善 和 退化 过 程 中 的 空间 分 布 


Fig.6 Spatial distributions of climate change, human activities and their combined 


effects on vegetation improvement and degradation 


(a) 气候 变化 


和 


LA 


行政 边界 
O 县 界 

O 乡镇 边界 
相对 贡献 率 /% 
mE 0-10 

EN 10-20 
EN 20-30 
E 30-40 
E 40-50 
EN 50-60 
EN 60-70 
EN 70-80 
EN 80-90 
EN 90-100 


FASS, 


0 10km 


(b) 人 类 活动 
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图 例 
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O 县 界 

O 乡镇 边界 
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EN 0-10 

EH 10-20 
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图 7 2000 一 2020 年 盐池 县 气候 变化 和 人 类 活动 在 植被 改善 区 的 相对 贡献 率 


Fig.7 Contribution rates of climate change and human activities to vegetation improvement in Yanchi County during 2000—2020 


漠 rcp e Ji HLA ER RE RAED 4 BPH 
KAI, ELAS APH BR 79035] LA SST oh 9H BG 16 E 
要 因素 ,尤其 在 典型 草原 分 布 区 域 , 达 到 91.67% 
( 表 3)。 


3.4 NPP 变 化 驱动 力 分 析 

3.4.1 植被 NPP 变 化 与 气候 因子 的 相关 性 气候 变 
化 是 引发 植被 NPP 变化 的 重要 原因 ,并 且 气 温和 降 
水 是 该 地 区 气候 变化 最 直接 和 敏感 的 因素 ”。 基 于 
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(a) 气候 变化 


0 10km 


图 8 2000—20204F. 


Fig. 8 Contribution rates of climate change and human activities to vegetation degradation in Yanchi County during 2000—2020 
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(b) 人 类 活动 


图 例 
行政 边界 

O BR 

O 乡镇 边界 
相对 贡献 率 /% 
EN 0-10 

Emu 10-20 
ES 20-30 
= 30-40 
= 40-50 
ES 50-60 
EU 60-70 
mmu 70-80 
mH 80-90 

0 10km EN 90-100 


大 水 坑 镇 


活动 在 植被 退化 区 的 相对 贡献 率 


表 2 气候 变化 与 人 类 活动 在 植被 不 同 变 化 程度 中 的 相对 贡献 
Tab.2 Relative contribution of climate change and human activities to different degrees of vegetation change 
a DR 气候 变化 人 类 活动 
eee eee eee 像 元 数 像 元 百分比 /9% 像 元 数 像 元 百分比 /% 
极 显著 改善 区 72008 68.91 17275 23.99 54733 76.01 
显著 改善 区 18004 17.23 14787 82.13 3217 17.87 
无 显著 改善 区 13844 13.25 8750 63.27 5058 36.73 
极 显 著 退 化 区 57 0.05 17 28.83 40 71.7 
显著 退化 区 41 0.04 10 24.39 31 75.61 
无 显著 退化 区 537 0.52 80 14.90 457 85.10 
表 3 气候 变化 与 人 类 活动 在 不 同 植被 类 型 改善 和 退化 中 的 相对 贡献 
Tab. 3 Relative contribution of climate change and human activities to the improvement 
and degradation of different vegetation types 
— 植被 改善 区 植被 退化 区 
气候 变化 相对 贡献 率 /% 人 类 活动 相对 贡献 率 /% 气候 变化 相对 贡献 率 /% 人 类 活动 相对 贡献 率 /% 

荒漠 32.02 67.98 29.47 70.53 
荒漠 草原 48.03 51.97 24.18 75.82 
典型 草原 69.61 30.39 8.33 91.67 
粮食 及 经 济 作物 64.83 35.17 28.57 71.43 
EX 61.37 38.63 - - 
落叶 灌 从 61.28 38.72 - = 


IR TER EXT 2000—2020 4E4E35] ANPP 与 年 降水 量 和 
了 相关 性 分 析 并 对 其 相关 性 系数 进行 


年 均 气 温 进 和 


显著 性 检验 ,将 结 


吉 果 划分 为 极 显著 正 / 负 相 关 (P< 


0.01) .显著 正 / 负 相 关 (0.01< P< 0.05) .无 显著 正 / 负 
相关 (P>0.05)6 个 级 差 区 域 ( 图 9)。 如 图 9 所 示 , 研 
究 区 植被 NPP 整体 上 与 降水 量 之 间 呈 正 相 关 。 其 
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(a) 气温 
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图 例 
行政 边界 

加 县 界 

L3 乡镇 边界 
ANPP 与 降水 量 的 
相关 性 

mm 极 显著 负 相关 


Em 极 显 著 正 相关 0 10km Em 极 显 著 正 相关 
9 2000 一 2020 年 盐池 县 植被 NPP 与 气候 因子 相关 系数 的 空间 分 布 
Fig.9 Spatial distributions of correlation coefficients between vegetation NPP and 
climate factors in Yanchi County during 2000—2020 
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中 , 星 显著 及 极 显著 正 相关 的 区 域 分 别 占 研 究 区 的 
63.37% 和 14.55%; 尽 管 研究 区 植被 NPP 与 气温 整体 
上 呈正 相关 (94.17% ) ,但 呈 显 著 及 极 显 著 正 相关 的 
区 域 分 别 仅 占 研究 区 的 7.85% 和 0.11%。 表 明 盐 池 
县 植被 对 降水 较为 敏感 ,降水 是 影响 盐池 县 植被 改 
善 的 主要 气象 因素 。 

3.4.2 人 类 活动 对 植被 NPP 的 影响 如 图 7、 图 8 所 
示 , 人 类 活动 主导 植被 改善 的 区 域 占 植 被 改善 区 面 
积 的 40.23%, 而 91.77% 的 植被 退化 区 则 由 人 类 活动 
所 主导 ,说 明 一 方面 人 类 活动 对 植被 改善 有 积极 作 
用 ,而 男 一 方面 表现 为 人 类 活动 驱动 了 局 地 的 植被 
退化 。 由 图 10a 可 以 统计 得 出 ,HNPP>0 的 面积 占 研 
究 区 总 面积 的 70.74% ,表明 人 类 活动 对 植被 的 负 疝 
影响 分 布 较为 普遍 ;HNPP<0 的 面积 仅 占 29.26% , E 
要 分 布 于 荒漠 及 荒漠 草原 地 区 .麻黄 山地 区 及 各 乡 
镇 的 灌溉 区 ,表明 人 类 活动 仅 在 这 些 区 域 对 植被 改 
善 起 促进 作用 。 

为 进一步 分 析 HNPP 的 变化 规律 ,将 HNPP 的 
正 负 与 HNPP 的 变化 趋势 进行 到 加 分 析 , 重 分 类 得 
到 4 种 结果 : 正 正 、 正 负 、 负 正 、 负 负 ( 图 10b)。 正 正 
表示 人 类 活动 对 植被 的 退化 作用 增强 , 占 总 面积 的 
9.31%; 正 负 表 示人 类 活动 对 植被 的 退化 作用 减弱 ， 
分 布 范围 最 广 (61.49% ); 负 正 表 示人 类 活动 对 植被 


的 改善 作用 减弱 ,分 布 范 围 最 小 (2.08%); 负 人 负 表 示 
人 类 活动 对 植被 的 改善 作用 增强 , 占 总 面积 的 
27.12%。 综 合 来 看 ,人 类 活动 对 研究 区 全 域 的 退化 
作用 都 有 所 减弱 ,并 在 亮 漠 及 荒漠 草原 地 区 表现 出 
显著 的 改善 作用 。 


4 讨论 


4.1 植被 时 空 变化 趋势 

近 20 a 来 ,研究 区 实施 退耕 还 林 还 草 及 封 山 禁 
牧 前 后 ANPP 的 变化 总 体 上 以 植被 改善 为 主 ,但 会 
因 和 气候 波动 出 现 小 幅 负 向 波动 ;植被 NPP 改 善 区 面 
积 比 例 远 大 于 NPP 退 化 区 面积 比例 。 本 研究 模拟 
植被 NPP 多 年 平均 值 分 布 范 围 为 14.75~514.07 g C- 
ma ,均值 为 202.00 g C-m”…a'。 这 与 倪 向 南 等 ”i、 
朱玉 果 等 ” 穆 少 杰 等 "、 李 柏 延 等 “利用 同类 方 
法 在 临近 区 域 的 估算 结果 接近 ,说 明 估 算 结 果 具 有 
较 好 的 置信 和 度 。 受 退耕 还 林 3、 封 山 禁 牧 玉 农业 
绿洲 扩张 引 \ 农 业 生 产 条 件 改善 等 因素 的 影响 ， 
人 类 活动 的 破坏 有 所 减少 ,大 量 坡 耕 地 被 转化 为 草 
地 和 林地 ,使 部 分 荒漠 草原 re UE BEI DX BAL 
ANPP 得 到 相应 提高 。 
4.2 植被 动态 变化 驱动 机 制 分 析 

降水 是 驱动 研究 区 植被 改善 的 主要 气候 因子 ， 
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图 10 2000 一 2020 年 盐池 县 年 均 HNPP 的 空间 分 布 及 人 为 活动 对 NPP 影 响 的 方向 与 趋势 
Fig. 10 Spatial distributions of annual average HNPP in Yanchi County during 2000—2020, 


and the direction and trend of influence of human activities on NPP 


且 植 被 NPP 与 降水 量 呈 显著 和 极 显 著 相 关 的 面积 
远大 于 其 与 气温 呈 显 著 和 极 显著 相关 的 面积 。 这 
与 周 伟 等 ” 的 关于 温 性 荒漠 草原 NPP 与 降水 量 的 
相关 性 最 显著 的 论断 一 致 ,主要 与 当地 近 30 a 的 降 
水 格局 转变 5 有关。 尽管 植被 恢复 对 黄土 高 原 局 
地 降水 量 的 增加 具有 积极 效应 ,但 大 规模 生态 建 
设 工 程 的 实施 ,使 得 区 域 植 物 生 态 需 水 量 增加 ,在 
降水 波动 背景 下 会 出 现 植 被 NPP 下 降 的 现象 ,这 印 
证 了 本 研究 中 降水 影响 植被 波动 变化 .促使 植被 稳 
定 增 长 的 结论 。 气 候 的 波动 会 对 植被 NPP 产 生 较 
大 影响 ”™, 本 研究 中 ,植被 NPP 随 气候 波动 呈 波 动 
中 上 升 趋势 。 已 有 研究 表明 ,本 研究 区 可 利用 年 降 
水 量 以 1.5 mm (10a) 的 降幅 减少 中 ,而 气温 呈 升 高 
趋势 ,这 势必 会 造成 植被 燕 散 量 增加 所 。 尽 管 气 
温 上 升 可 延长 植物 生长 季 , 使 植被 NPP 增 加 3, 但 
高 温 会 加 剧 植被 蒸发 蒸腾 ,降低 水 分 利用 效率 ,其 
至 出 现 土壤 干 层 , 限 制 了 草地 植被 生长 的 。 本 人 研究 
中 ,气温 与 植被 ANPP 总 体 上 呈正 相关 ,但 相关 性 不 
显著 ,说 明 植被 改善 的 减缓 可 能 会 受到 气温 升 高 的 
抑制 作用 。 

研究 区 退耕 还 林 、 退 牧 还 草 工程 的 实施 及 草原 
禁 牧 封 育 政策 的 落实 ,对 土壤 侵蚀 控制 .植被 覆盖 


状况 改善 均 起 到 了 积极 作用 '* ,表现 为 正 向 人 类 干 
扰 增 加 而 负 向 人 类 干扰 减少 ,使 荒漠 草地 退化 现象 
得 到 进一步 缓解 ,植被 状况 得 到 显著 改善 。 但 在 空 
间 分 布 上 具有 一 定 的 空间 异 质 性 ,具体 表现 为 : 
40.23% 的 植被 改善 区 由 人 类 活动 为 主 所 驱动 ,而 
91.77% 的 植被 退化 区 则 由 人 类 活动 为 主 所 驱动 。 
但 由 于 退耕 还 林 还 草 . 退 牧 还 草 .于 封禁 牧 等 生态 
恢复 政策 的 实施 ,人 类 活动 对 植被 退化 区 NPP 的 
负面 影响 逐渐 减弱 。 
43 气候 变化 与 人 类 活动 对 植被 变化 的 相对 贡 
通过 NPP 差 值 比较 法 ,分 离 气候 变化 与 人 类 活 
动 对 NPP 的 影响 , 既 可 以 定量 计算 NPP 的 人 为 影响 
值 ,评价 生态 工程 的 实施 成 效 ,也 可 以 厘定 气候 变 
化 对 NPP 影 响 的 贡献 率 。 本 研究 发 现 气候 变化 和 
人 类 活动 及 其 二 者 共同 作用 存在 显著 的 空间 异 质 
PA Te Ep Ply a As E A] H He SEER a i TT 
域 植被 变化 的 归 因 结果 一 致 。 本 研究 中 ,负面 的 人 
类 活动 如 城市 扩张 .矿产 资源 开采 、 偷 牧 等 依旧 是 
研究 区 植被 退化 的 主要 驱动 力 , 但 范围 较 小 ,而 围 
栏 放牧 、 禁 牧 .退耕 还 林 和 生态 修复 工程 等 对 大 范 
围 的 植被 改善 有 积极 作用 。 


202208.00050v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


FHEWY 


5 结论 


本 研究 利用 2000—2020 ^F. MODIS NDVI fit [H] 

列 数据 、 气 象 数据 和 地 表 和 覆盖 数据 ,以 NPP VE 
植被 变化 的 衡量 指标 ,基于 像 元 尺度 通过 CASA 模 
型 和 Thornthwaite 纪 念 模型 模拟 盐池 县 因 生 态 工程 
的 实施 对 区 域 NPP 时 空 变化 的 影响 ,并 通过 NPP 25 
i 1 aie Ae re eh 
植被 恢复 和 退化 的 相对 贡献 率 。 得 出 以 下 结 
(1) 2000 一 2020 年 盐池 县 植被 ANPP aias " 
增加 趋势 。 气 候 变化 和 人 类 活动 共同 驱动 了 盐池 
县 植被 NPP 的 总 体 增加 趋势 (2) 气候 变化 和 人 类 
活动 在 植被 改善 和 退化 中 的 相对 贡献 率 表现 出 显 
著 的 空间 异 质 性 。 其 中 ,在 植被 改善 区 ,气候 变化 
的 相对 贡献 率 为 59.77% ,而 人 类 活动 的 相对 贡献 率 
为 40.23%; 在 植被 退化 区 ,人 类 活动 的 相对 贡献 率 
为 91.77% ,远大 于 气候 变化 的 相对 贡献 率 。 说 明 气 
候 变化 和 人 类 活动 一 方面 促进 盐池 县 总 体 上 的 植 
被 恢复 ,同时 也 造成 极 小 部 分 区 域 植被 退化 ,但 退 
化 趋势 显著 减弱 。 
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Effects of climate change and human activities on 
net primary productivity of vegetation in Yanchi County 


LI Shiyao', CONG Shixiang, WANG Rongrong, YU Hailong, HUANG Juying 
(1. School of Geography and Planning, Ningxia University, Yinchuan 750021, Ningxia, China; 2. School of Ecology and 


Environment, Ningxia University, Yinchuan 750021, Ningxia, China) 


Abstract: This study aims to quantitatively differentiate the relative contribution rate of climate change and 
human activities during the vegetation cover change in the typical agricultural pastoral ecotone of Northern China 
under the background of the "grain to green project". Taking vegetation change from 2000 to 2020 in the Yanchi 
County, Ningxia Province as the study object and using MODIS13Q1-NDVI, land cover data, and meteorological 
data from 2000 to 2020, the annual potential net primary productivity and actual net primary productivity were 
simulated using the models of Thornthwaite Memorial and Carnegie Ames Stanford Approach. The spatial and 
temporal change characteristics of the net primary productivity (NPP) in the Yanchi County and its driving force 
were analyzed, and the relative contribution rates of meteorological factors and human activities in the process of 
vegetation cover change were quantitatively calculated using trend analysis, correlation analysis, and difference 
comparison. The results showed the following: (1) The NPP totally presented an increasing trend, but there was 
an obvious spatial heterogeneity in the NPP change trend and degree, mainly manifested when the vegetation 
NPPs of the improved areas were far greater than the those of the degraded areas, and the extents of improvement 
and degradation showed significant spatial heterogeneity. The vegetation-improved areas were mainly distributed 
in the desert and desert steppe (where the "grain to green project" was conducted) and in irrigated areas, whereas 
the vegetation- degraded areas were mainly distributed in the edge zones of the desert and desert steppe. (2) The 
vegetation change attribution analysis results indicated that in the vegetation- improved areas, the relative contribution 
rate of climate change (59.77%) was higher than that of human activities (40.23%). Meanwhile, in the vegetation- 
degraded areas, the relative contribution rate of human activities (91.7796) was higher than that of climate change 
(8.23%). (3) The driving force analysis results indicated that the vegetation NPP was positively correlated with 
precipitation but weakly correlated with temperature. Human activities were the main driving force that drove the 
NPP decrease in the vegetation- degraded areas. In summary, climate change was the main driving force that drove 
the NPP improvement in the vegetation- improved areas, whereas human activities were the main driving factors 
that led to the NPP decrease in the vegetation- degraded areas. Climate and human activities jointly dominated the 
improvement of the total ecological environment in the Yanchi County. 


Key words: net primary productivity (NPP); desert steppe; vegetation change; driving force; Yanchi County 


